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Abstract 

Quintessence issues can be achieved by taking into account higher order curva- 
ture invariants into the effective action of gravitational field. Such an approach is 
naturally related to fundamental theories of quantum gravity which predict higher 
order terms in loop expansion of quantum fields in curved space-times. In this 
framework, we obtain a class of cosmological solutions which are fitted against 
cosmological data. We reproduce encouraging results able to fit high redshift su- 
pernovae and WMAP observations. The age of the universe and other cosmological 
parameters are discussed in this context. 
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1 Introduction 



The existence of a dark energy term into cosmological dynamics has become a paradigm 
in the last five years since several observational campaigns give reliable indications on 
the apparent acceleration of the universe. 

High redshift supernovae surveys [1, 2], CMBR data [3, 4, 5], Sunyaev-Zeldovich / X-ray 
methods [6] and other approaches provide a new picture of the universe. It can be repre- 
sented as an isotropic, homogeneous, spatially flat 4-dim manifold filled with about 30% 
of baryonic and non-baryonic matter and about 70% of dark energy, simply referred as 
"cosmological component" . 

Also the very recent WMAP data [7] confirm such a picture with extremely low uncer- 
tainties in the estimate of cosmological parameters. 

It is evident that the cosmological component should be the ingredient capable of gener- 
ating the accelerated expansion, but, till now, its real nature is a puzzle which seems far 
to be solved. 

Many approaches has been developed. Cosmological constant is the most straightforward 
candidate, however it is ruled out since its observational value (constrained by observa- 
tions) differs of 120 order of magnitude from the theoretical prediction of QCD [8, 9, 10]. 
This argument does not allow to interpret cosmological constant as the vacuum energy 
of gravitational field unless an evolutionary mechanism is invoked in order to explain 
dynamics starting from the huge early values of energy {cosmological constant problem). 
Besides, the observed comparable amounts (in order of magnitude) of matter and dark 
energy sets a strong fine-tuning problem which cannot be overcome considering wide 
ranges of initial data {coincidence problem) [11]. 

A second approach is to consider the cosmological component as a dynamical term. This 
scheme, usually called quintessence, can be achieved adding scalar fields into Einstein 
gravity. Such scalar fields are assumed rolling their interaction potential [11, 12], as 
in the case of inflation but in a very different energetic regime (today instead of early 
universe). However major shortcomings come out due, essentially, to the ad hoc forms 
of self-interaction potential. Several forms of potential (inverse power law, exponential, 
etc, [9, 11]) achieve quintessence prescriptions (e.g. accelerated behaviour, Q\ ~ 0.7, 
coincidence problem, etc.) but none of them seems to be directly related to some funda- 
mental quantum fleld theory. The situation is, somehow, similar to that of inflationary 
cosmology: inflation is a good paradigm but no single model matches all requests. 
On the other hand, alternative approaches can be pursued starting from some funda- 
mental theory [13, 14, 15]. These schemes aim to improve the quintessence approach 
overcoming the problem of scalar fleld potential, generating a dynamical source for dark 
energy as an intrinsic feature. The goal would be to obtain a comprehensive model ca- 
pable of linking the picture of early universe to the today observed one; that is, a model 
derived from some effective theory of quantum gravity which, through an inflationary 
period results in the today accelerated Priedmann expansion driven by some f^A-term. 
Prom this point of view, theories including torsion [16, 17, 18] or higher order curvature 
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invariants [19] naturally come into the game. In fact, the former ones allow to include 
spin matter fields at a fundamental level in General Relativity, while the latter ones come 
out in every quantization scheme of matter fields in curved space-time [9, 20, 21]. Specif- 
ically higher-order terms in curvature invariants as R^, R^j^tjR'^" , RDR are inescapable if 
we want to obtain an effective action of gravity closed to the Planck epoch. This scheme 
naturally give rise to inflationary behaviours [22, 23, 24]. 

In [19], it is investigated the possibility that such terms could act also today as a non- 
clustered form of dark-energy, providing accelerated behaviours. Besides, other authors 
have recently discussed curvature self-interactions of cosmic fluid to get cosmological con- 
stant [25]. 

In this paper, we want better to detail such an approach. We provide exact solutions and, 
following the scheme in [18], we try to obtain observational constraints on the so called 
Curvature Quintessence in order to see if curvature contributions can actually match the 
data of recent surveys. 

The paper is organized as follows. In Sec. 2, we summarize the curvature quintessence 
approach and derive a class of exact cosmological solutions. Sect. 3 is devoted to the 
matching with observational data. In particular, we fit our solutions against the SNIa 
data and derive the age of the universe. A further discussion is carried out in Sec. 4 
considering the very recent WMAP data, which seem to better constrain the parameters 
of cosmological solutions. Sec.5 is devoted to the conclusions. 



2 Curvature Quintessence 

A generic fourth-order theory of gravity, in four dimensions, is given by the action [19], 

A= d'^X^/^ \f{R) + jC{matter)\ , (1) 



where f{R) is a function of Ricci scalar R and C(matter) is the standard matter Lagrangian 
density. We are using physical units SttGn — c — h — 1. The field equations are 

which can be recast in the more expressive form 

n — f? ^ „ TD — T^(c«''^') I rp{rnatter) /o\ 

where an stress-energy tensor has been defined for the curvature contributes 

T^7'^ = [la., [fiR) - Rf{R)] + fiRy^'iga^g,. - ga^g,.)] (4) 
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and 



rp{matter) 



^ rj^{matter) 



(5) 



is the stress-energy tensor of matter. We have taken into account the nontrivial couphng 
to geometry; prime means the derivative with respect to R. If /(-R) = it! + 2 A, we recover 
the standard second-order gravity. 

In a Priedmann- Robertson- Walker (FRW) metric, the action (1) reduces to the point-hke 
one: 



(curv) 



dt 



C{a,a;R,R) + C^rnatter) 



(6) 



where the dot means the derivative with respect to the cosmic time. In this case the 

scale factor a and the Ricci scalar R are the canonical variables [26, 27]. 

It has to be stressed that the definition of R in terms of a, d, a introduces a constraint 

in the action (6) [19], by which we obtain by the Lagrange multiplier technique the 

lagrangian 



>C = C^curv) + Ci^matter) = [f{R) - Rf{R)] + Qaa^f'{R)+ 

+6a^dRf"{R) - 6kaf'{R) + a^PirnaUer) , 



(7) 



(the standard fluid matter contribution acts essentially as a pressure term [28]). The 
Euler-Lagrange equations coming from (7) give the second order system: 



21"- 



k 



-Pitot) 



and 



f{R) |i? + 6 
constrained by the energy condition 



a fd\ k 



a 



aj 



a" 



0, 



ay k _1 

Using Eq.(lO), it is possible to write down Eq.(8) as 

'd\ 1 
.a/ 6 



P{tot) + ^P(tot) 



(9) 



(10) 



(11) 



The accelerated or decelerated behaviour of the scale factor depends on the r.h.s. of (11). 
The accelerated behaviour is achieved if 

P(tot) + Sp{tot) < . (12) 

To understand the actual effect of these terms, we can distinguish between the matter 
and the geometrical contributions 



P{tot) — P{curv) + P{matter) P{tot) — P(curv) + P(matter) ■ 



(13) 
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Assuming that all matter components have non-negative pressure, Eq.(12) becomes: 



P{curv) > ^Pifot) ■ (14) 

The curvature contributions can be specified by considering the stress-energy tensor (4); 
we obtain a curvature pressure 

Picurv) = 77^ {2 (^) Rf'iR) + Rf'iR) + R'f'iR) - \ m) - Rf'iR)]] , (is) 

and a curvature energy-density: 

P(curv) = [/(^) - ^/'(^)] - 3 ^/"(^)} ' (16) 

which account for the geometrical contributions into the thermodynamical variables. 
It is clear that the form of f{R) plays an essential role for this model. We choose the 
f{R) function as a generic power law of the scalar curvature and we ask for power law 
solutions of the scale factor. 
Summarizing, we consider: 

f(R)^foR\ a(^) = ao(^)^ (17) 

The interesting cases are for a > 1 which give rise to accelerated expansion. 
Let us now concentrate now concentrate on the case with P(matter) = 0. Inserting Eqs.(17) 
into the above dynamical system, for a spatially fiat space-time {k = 0), we obtain an 
algebraic system for the parameters n and a 

' a[a{n - 2) + 2/1^ - 3n + 1] = 

(18) 

_ a{n^ + a{n - 2 - n - \)\ = n{n - l)(2n - 1) 



from which the allowed solutions: 

a = n = 0, 1/2, 1 

2n2-3n + 1 , ,^ ^^^^ 

a — , Vn but n 2 . 

2 — n 



The solutions with a = are not interesting since they provide static cosmologies 
with a non evolving scale factor^. On the other hand, the cases with generic a and n 
furnish an entire family of significative cosmological models. By the plot in Fig.l we see 
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Figure 1: The behaviour of a in term of n. It is evident a region in which the power of the 
scale factor is more than one and a region in which it is always negative. The plot on the right 
represent the values of n between and 1.5. 



that such a family of solutions admit negative and positive values of a which give rise to 
accelerated behaviours. 

Using Eqs.(15) and (16) we can also deduce the state equation (the barotropic index) 

for the family of solution a — — We have 

2 — n 



/6n2-7n-l\ 
= - -9n + 3 j ' (2°) 

which clearly is W(^curv) — 1 forn — > oo. This fact shows that the approach is compatible 
with the recovering of a cosmological constant. 

Eq.(20) is plotted in Fig.2. From the two plots, we can observe that the accelerated 
behaviour is allowed only for W(^curv) < as requested for a cosmological fluid with 
negative pressure. 

The whole approach seems intriguing. In fact we are able to describe the accelerated 
phase of universe expansion simply as an effect of higher order curvature terms which 
provide an effective negative pressure contribute. In order to see if such behaviour is 
possible for today epoch we have to match our model with observational data. 

^This result match with the standard General Relativity case (n=l) in absence of matter. 
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Figure 2: Behaviour of Wcurv against n. We have drawn a Une for the value Wc 
emphasizing the cosmological constant state equation value. 
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3 Matching with SNIa observations 

To verify if the curvature quintessence approach is an interesting perspective, we have 
to match the model with the observational data. In this way, we can constrain the 
parameters of the theory to significant values. First we compare our theoretical setting 
with the SNela results. As a further analysis, we check also the capabihty of our model 
in the universe age prediction. 

We have stressed in Sec.l that SNela observations have represented a cornerstone in the 
recent cosmology, pointing out that we live in an expanding accelerating universe. This 
result has been possible in relation to the feature of supernovae to be considered standard 
candles via the Phillips amplitude-luminosity relation. 

To test our cosmological model, we consider supernovae observations reported in [1] [2]. 
Wc have compiled a combined sample of 79 SNela discarding 6 likely outliers SNela as 

discussed in [1]. 

Starting from these data, it is possible to perform a comparison between the theoretical 
expression of the distance modulus 

/x(^) = 51og^dz,(z) + 25, (21) 
-no 

and its experimental value for SNela, z is the redshift (c is the light speed, hereafter we 
will use standard units). In general, the luminosity distance dL can be expressed as: 
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where E{() = —. From (17), the Hubble parameter is 
Ho 



H{z)^Ho{l + z)'^\ (23) 



where a depends on n as in Eq.(19). A fit by SNela data provides significant values of 
the parameter n which assigns the curvature function. 

The best fit is performed as in [29] minimizing the calculated between the theoretical 
and the observational value of distance modulus: 

X [Hq, n) - 2^ (24) 

i IJo,i ' mz,i 



for the above considered data [1, 2]. 
Now, the luminosity distance becomes: 

Hq Jo 

which, after integration, gives 



(25) 



dLiz,Ho,n) = 

riQ \cx — 1 



1 + ^) 



'l^z)^^ - 1 



(26) 



This expression is not defined for ck = 0, 1 which physically correspond to static universes 
and Milne ones; such models are not interesting for our purposes. The range of n can be 

divided into intervals taking into account the existence of singularities in (26). Thus, the 

11 1 
fit is performed in five intervals of n, which are: n < -{1 — -\/3) , -(1 — -\/3) < n < - , 

^ < n < 1 , 1< n< ^(1 + V3) , n > ^(1 + \/3). 

In order to define a limit for Ho^ we have to note that the Hubble parameter, as a function 
of n, has the same trend of a. showed in Fig.l. We find that for n lower than -100, the 
trend is strictly increasing while for n positive, greater than 100, it is strictly decreasing. 
In relation to this feature, we have tested the values of n ranging in these limits because, 
as we shall see below, out of this range, the value of the age of universe becomes mani- 
festly non-physically significant. The results of the fit are showed in Table 1. In Fig.3, 
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Range 








-100 <n< 1/2(1 - V3) 


65 


-0.73 


1.003 


1/2(1-^3) < n < 1/2 


63 


-0.36 


1.160 


1/2 < n< 1 


100 


0.78 


348.97 


l<n< 1/2(1 + ^3) 


62 


1.36 


1.182 


1/2(1 + V3) < n< 3 


65 


1.45 


1.003 


3 < n< 100 


70 


100 


1.418 



Table 1: Results obtained by fitting the curvature quintessence models against SNela data. 
First column indicates the range of n, column two gives the relative best fit value of Hq, column 
three n^"^*'*, column four the index. 



we present the different contour plots for each evaluated range. 



From Tab.l and Fig.3, it is evident that only in the ranges 1/2(1 — ^/3) < n < 1/2 and 
l/2<n<lwe can achieve a constraint on the values of n. Besides, we can exclude the 
range l/2<n<las physically not interesting, in relation to the best fit value of Hq 
and the result. In the other cases, the test gives significant best-fit values both for 
Hubble parameter and x^. But the contour plots are completely degenerate with respect 
to the n parameter. This occurrence suggests that varies slightly with respect to n in 
this cases, hindering the possibility to constrain n. 

In order to discriminate the significant value of n, we have also performed a test of the 
model in relation to the capability of estimating the age of the universe. 

The age of the universe can be simply obtained, from a theoretical point of view, if 
one knows the value of the Hubble parameter. In our case, from the definition of if, and 
using the relations (17), we have: 

t - aH-^ (27) 



which is 



The age of universe can be obtained substituting H — Hq in (28). 
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Figure 3: Contour plots for the considered n-ranges using the SNIa data. It is evident that in 
the interval 1/2 < n < 1 the models are non physical. The ranges 1/2(1 - V3) < n < 1/2 and 
1/2 < n < 1 give some indications on the best fit value of n, in the other cases the matching 
with SNIa is completely degenerate. 9 



Range 


AH{kms-^Mpc-^) 


An 


t{n^^'*){Gyr) 


-100 <n< 1/2(1 - V3) 


50-80 


-0.67 <n< -0.37 


23.4 


1/2(1 - V3) < n < 1/2 


57-69 


-0.37 <n< -0.07 


15.6 


l<n< 1/2(1 + ^3) 


56-70 


1.28 <n< 1.36 


15.3 


1/2(1 + ^3) <n<2 


54-78 


1.37 <n< 1.43 


24.6 



Table 2: The results of the age test. In the first column is presented the tested range. Second 

column shows the 3o"-range for Hq obtained by supernovae test, while in the third wc give the n 
intervals, i.e. the values of n which allow to obtain ages of the universe ranging between lOGyr 
and 18 Gyr. In the last column, the best fit age values of each interval are reported. 

We evaluate the age taking into account the intervals of n and the 3cr-range of variability 
of the Hubble parameter deduced by the Supernovae fit. We have considered, as good 
predictions, age estimates included between lOGyr and 18G'yr. 



By this test, we are able of refine the allowed values of n. The results are shown in Table 
2. First of all, we discard the intervals of n which give negative values of t. Eq. (27) 
shows that negative values of the age of the universe are obtained for negative values of 
a, so we have to exclude the range 1/2 < n < 1 and n > 2 (Fig.l). Conversely, the other 
ranges, tested by SNIa fit (Tab.l), become narrower, strongly constraining n. 



A further check for the allowed values of n is to verify if the interesting ranges of n 
provide also accelerated expansion rates. This test can be easily performed considering 

, using the relation (17) and 



the definition of the deceleration parameter qo — — ( t^t ) 



the definition of a in term of n. To obtain an accelerated expanding behaviour, the scale 
factor a{t) = aot" has to get negative or positive values of a greater than one. Wc obtain 
is that only the intervals —0.67 < n < 0.37 and 1.37 < n < 1.43 provide a negative 
deceleration parameter with a > 1. Conversely the other two intervals of Tab. 2 do not 
give interesting cosmological dynamics, being > and < ck < 1. 



4 WMAP Age test 

A further test of the model can be performed by the age estimate obtained by the 
WMAP campaign. The WMAP ( Wilkinson Microwave Anisotropy Probe) mission aims 
to determine the geometry, the content and the evolution of the universe through a full sky 
map of the temperature anisotropies of cosmic microwave background radiation [7] . The 
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first year observational results has been published since few months. These observations 
indicate as the best fit model a A-tcrm cosmological model with about 70%-content of 
cosmological origin, Hubble constant value 7lio o| km s~^Mpc~^ and an age estimate of 
13.7t',iGyr. 

Using this last result, we can improve the constraints on n in relation to the very low 
error (1%) of WMAP age estimator. 

We use the same approach of the previous section, the only difference is to consider as 
physically interesting only the age prediction ranging between 13.5Gyr and 13.9Gyr. The 
results are shown in Tab. 3. 

It is evident that this test narrows the range of physical interest. If we take into account 
also the capability of providing acceleration, we obtain that the interesting values of n 
are -0.450 <n< -0.370 and 1.366 <n< 1.376. 

These results could represent a selection for the allowed form of fourth-order gravity 
action as a power law of curvature Ricci scalar. 



Range 


AH{kms-'^Mpc-^) 


An 


qo 


-100 <n< 1/2(1 - 73) 


50-80 


-0.450 < n < -0.370 


< 


1/2(1-^3) <n< 1/2 


57-69 


-0.345 < n < -0.225 


> 


l<n< 1/2(1 + ^3) 


56-70 


1.330 <n< 1.360 


> 


1/2(1 + V3) < n< 2 


54-78 


1.366 <n< 1.376 


< 



Table 3: The results of the age test for the curvature quintessence model based on the prediction 
of WMAP observations. In the last cohimn wc show the sign of the deceleration parameter as 
a testify of the potential accelerating rate of the model. 



5 Conclusions 

In this paper we have analyzed a geometrical approach to quintessence given by con- 
sidering a fourth-order theory of gravity and we have matched the cosmological models 
derived with observational data. This scheme, proposed as curvature quintessence [19], 
represents an approach to describe an accelerated expanding universe dominated by a 
cosmological component without using scalar fields. 

We stress that such an approach has a natural background in several attempts of quantize 
gravity, because higher-order curvature invariants are generated to renormalize quantum 
field theories on curved space times [21]. 

Wc have studied curvature quintessence neglecting matter as a first significative approxi- 
mation. This approximation is possible if we consider that, at one loop level, matter and 
gravity have the same "weight" in early universe [20] and today the amount of cosmo- 
logical and matter components, at cosmological density level, are comparable. In both 
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regimes, the overall evolution can be assigned only by the cosmological component which 
gives a good approximation of dynamics. We ask for a power law functions for the action 
(in terms of Ricci curvature scalar) and for the scale factor (in terms of cosmological 
time). With these choices, we obtain a family of exact solutions so that the model is 
characterized by one parameter (specifically n). To check these solutions, we have fitted 
the model with observational data. A straightforward test is a comparison with SNIa 
observations referring to the well known data of SCP {Supernovae Cosmology Project [1]) 
and of HZT {High-Z search Team [2]). 

The model fits these data and provides a constrain on the parameter n. Unfortunately 
this test is not decisive to feature curvature quintessence since the likehhood curves are 
degenerate with respect to n. 

To improve this result, we have performed a test with the age of the universe. 
In a first case, our test has been conceived considering, as good estimates for age of uni- 
verse, values ranging between lOGyr and 18 Gyr. Considering that we take into account 
only accelerated dynamics, we deduced that for —0.67 <n< 0.37 and 1.37 < n < 1.43 
curvature quintessence is a model capable of mimicking the actual universe. 
In order to better refine these ranges, we have then considered a test based on WMAP 
age evaluation. 

In this case, the age ranges between 13.5Gyr and 13.9Gyr. This fact reduces the allowed 
intervals of the parameter n, which are —0.450 <n < —0.370 and 1.366 <n < 1.376. In 
conclusion, we can say that a fourth order theory of gravity of the form 

f{R) = foR'^' (29) 

with e ^ —0.6 or e ^ 0.4 can give rise to reliable cosmological models which well fit 

SNela and WMAP data. 

In this sense, we need only "small" corrections to second order Einstein gravity action in 
order to achieve quintessence issues. 

Indications in this sense can be found also in a detailed analysis of f{R) cosmological 
models performed against CMBR constraints, as shown in [30]. 
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